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Qh! A search for exotic unstable neutral and charged heavy leptons as well as for 

stable charged heavy leptons is performed with the L3 detector at LEP. Sequential, 
vector and mirror natures of heavy leptons are considered. No evidence for their 
existence is found and lower limits on their masses are set. 
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Introduction 

Many theories beyond the Standard Model of the electro-weak interation pP predict the exis- 
tence of new heavy leptons [2j with observable production cross sections and cleanly identifiable 
final states. Results on this subject obtained at the Z resonance by LEP and SLC experiments 
are found in Reference El Results obtained by the LEP experiments at centre-of-mass energies, 
y/s = 133—189 GeV are found in References HI and El Here we report on a search for exotic 
pair-produced heavy leptons. The data used in this analysis were collected with the L3 detec- 
tor [6 at LEP during 1999 at <Js = 192-202 GeV and 2000 at <Js = 200-208 GeV with a 
total integrated luminosity 450 pb _1 . The results are combined with our earlier data recorded 
at y/s = 133 - 189 GeV. 

The exotic leptons discussed in this paper may be classified by their SU{2) x U(l) quantum 
numbers |2j: 

• Sequential leptons: they exist in the simplest extension of the Standard Model where a 
fourth family with the same quantum numbers is added to the known fermionic spectrum. 
The indirect limits on existence of these particles can be found in Reference El 

• Vector leptons: these particles occur in both left-handed and right-handed weak isodou- 
blets. Among others, the E 6 group model jH] predicts them. 

• Mirror leptons [§]: these particles have chiral properties which are opposite to those 
of ordinary leptons, i.e. the right-handed components form a weak isodoublet and the 
left-handed ones are weak isosinglets. Mirror fermions appear in many extensions of the 
Standard Model and provide a possible way to restore left-right symmetry at the scale of 
the electroweak symmetry breaking. 

Production and decay of exotic leptons 

Charged heavy leptons, L*, are pair-produced through the s-channel via 7 and Z boson ex- 
change, while for heavy neutral leptons, L°, only the Z boson exchange is present. The total 
cross sections are in the range of 1-3 pb at masses well below the beam energy, E^ eam ^ and fall 
as the mass of the lepton approaches the beam energy. The neutral heavy leptons are either 
of Dirac or Majorana types. The main difference between pair-produced Dirac and Majorana 
neutral heavy leptons is the dependence of the cross section on their velocity /3. A (3(3 — (3 2 ) 
term is present in the Dirac case while the cross section in the Majorana case is proportional to 
/3 3 . This implies that the Majorana cross section falls more rapidly with mass than the Dirac 
one. In this search, heavy leptons couple to the electron, muon, or tau families. Only the case 
where both heavy leptons decay into the same light lepton family is considered. 

Exotic leptons decay via the charged or neutral weak currents through the mixing with a 
light lepton, except for sequential leptons, which decay only through charged current: 

I± _> Ui w ± or L ± -> £±Z 

L° -► £ ± W T or L° -> v{L 

where £ = e, //, r. For masses, mi,, below 100 GeV, heavy leptons decay predominantly through 
the W boson. The branching fraction of L° — > eW is about 90% for m^ = 100 GeV and 100% 
for rriL < 90 GeV. Therefore, only the decay into the W boson is considered. 



The decay amplitude of L° — > ^W^ contains a mixing parameter U . The mean decay 
length, D, as a function of \U\ 2 and mass is given by D = fi^cr^ oc l3\U\" 2 mf j [TT], where tl 
is the lifetime of the heavy lepton and a ~ —6. This implies the decay can occur far from 
the interaction point if the particle has a low mass or a very small coupling. To ensure high 
detection and reconstruction efficiencies, this search is restricted to decay lengths below 1cm, 
which limits the sensitivity to the mixing parameter to \U\ 2 > O(10~ n ). 

Two additional possibilities for the charged heavy lepton decays are considered: 

• The charged leptons decay through a weak charged current interaction that conserves 
lepton number, L+ — >L°W+. 

• The charged leptons are stable. This is the case if the associated neutral lepton is heavier 
than its charged partner and the coupling with light leptons is negligible. 

Event simulation 

The generation of heavy leptons and their decay is performed with the TIPTOP ^2] Monte 
Carlo program which takes into account initial state radiation and spin effects. A mass range 
between 50 and 100 GeV for the heavy leptons is considered. For the simulation of back- 
ground from Standard Model processes the following Monte Carlo programs are used: KK2f J3] 
(e + e~ -» qq(7)), PYTHIA HI] (Ze+e", ZZ), BHWIDE HSj (e+e" -> e+e"^)), KORALZ [IE] 
(e+e" -> ;U + /!~(7) and e+e" -> r+i-fr)), KORALW [H] (e+e" -> W+W"), PHOJET [IB] 
(e+e~ -> e+e-qq), DIAG36 [E] (e+e" -> e + e~£ + £-), and EXCALIBUR f2UJ for other four- 
fermion final states. The Monte Carlo events are simulated in the L3 detector using the 
GEANT3 program [21], which takes into account the effects of energy loss, multiple scattering 
and showering in the materials. Time dependent detector inefficiencies, as monitored during 
the data taking period, are also simulated. 

Search for unstable heavy leptons 

Lepton and jet identification 

Unstable heavy lepton signatures are characterized by hadronic jets, isolated leptons, and 
missing energy and transverse momentum. 

An electron is identified as a cluster in the electromagnetic calorimeter with an energy larger 
than 4 GeV matched to a track in the (R, 0) plane to within 20 mrad. The cluster profile should 
be consistent with that expected for an electron. The polar angle of the electron candidate must 
satisfy | cos#| < 0.94. A muon track is required to have track segments in at least two out of 
three layers of muon chambers and to point back to the interaction region. The muon should 
have a momentum greater than 4 GeV and a polar angle defined by | cos#| < 0.8. The lepton 
is isolated if less than 5 GeV is deposited in a 30° cone around its direction. Moreover, to 
eliminate photon conversions, only one matched track is required for isolated electrons. Jets 
are reconstructed from electromagnetic and hadronic calorimeter clusters using the Durham 
algorithm (22| with a jet resolution parameter y cut = 0.008. The jet momenta are defined by 
the vectorial energy sum of calorimetric clusters. 

The most important backgrounds are W+W~, qq(7) and ZZ production. The W+W~ 
process is rejected by cuts on visible mass, number of hadronic jets, jet angles and energies 



and lepton energies. The qq(7) process is suppressed by demanding transverse momentum 
imbalance, no energy in the very forward-backward regions, and a missing momentum vector 
pointing into the detector. For ZZ event rejection, a cut on the energy sum of the two isolated 
leptons is applied. 

Search for unstable neutral heavy leptons 

The search for pair-produced neutral heavy leptons requires two isolated leptons (e, //, or 
r) of the same family together with the decay products of W bosons, i.e. e + e~ ^L°L — > 
£ + £~W + W~ . Hadronic events with visible energy greater than 60 GeV and charged track 
multiplicity greater than 3 are considered. 

For the case where both neutral heavy leptons decay to either electrons or muons, events 
must also satisfy the following criteria: 

• The number of reconstructed jets plus isolated leptons is at least 3. 

• The event contains at least two isolated electrons or two isolated muons. Figure 1 shows 
the energy in a 30° cone for the second most energetic electron candidate. 

• The energy sum of the two isolated electrons or muons must be less than 0.7x Eb eam . 
Figure 2 shows the sum of the energies of two muons. 

After this selection, 11 events remain in the data for the electron decay mode while 10.7 ± 0.5 
background events are expected. For the muon decay mode, 5 events remain in the data while 
2.3 ± 0.2 background events are expected. 

For the case where both neutral heavy leptons decay to tau leptons, each of the taus decays 
independently into hadrons, muons, or electrons. When both tau leptons decay into either 
muons or electrons, the above selection is applied with the exception that we allow the isolated 
leptons to be either two electrons, two muons, or one muon and one electron. For the case in 
which at least one of the tau decays into one charged hadron, events are required to satisfy the 
the following criteria: 

• The number of reconstructed jets plus isolated leptons is at least 4. 

• The polar angle of the missing momentum must be in the range 25° < miss < 155°, 
as depicted in Figure 3, and the fraction of the visible energy in the forward-backward 
region, 9 < 20° and 9 > 160°, must be less than 40%. 

• The angle between the most isolated track and the track nearest to it must be greater 
than 50°, or the angle between the second most isolated track and the track nearest to 
it must be greater than 25°. The transverse momenta of the two most isolated tracks 
must be greater than 1.2 GeV, and at least one track must have a transverse momentum 
greater than 2.5 GeV. 

• The electron and muon energies must be less than 50 GeV. 

After applying the above selection and also accounting for the case where both taus decay into 
electrons or muons, 114 events remain in the data, while 121.5 ±1.6 background events are 
expected. 



Search for unstable charged heavy leptons 

Decay into a light neutrino, L 1 * 1 — ► z/^W ± 

For this search, two sets of cuts are used to identify the topology with one hadronic and one 
leptonic W-decay as well as two hadronic W-decays. For both selections, the visible energy is 
required to be greater than 50 GeV, the multiplicity of charged tracks to be greater than 3, and 
the fraction of the total visible energy in the forward-backward region to be less than 25%. 

For the decay L + L~ — > v^pn W + W~ — ► i/£Ug£i/gqq' , events satisfying the following criteria 
are selected: 

• The event contains at least one isolated electron or muon with energy greater than 15 GeV 
and less than 60 GeV. 

• The number of reconstructed jets plus isolated leptons is at least 3. 

• The polar angle of the missing momentum must be in the range 25° < 9 m i SS < 155°. 

• The sum of the energies of the hadronic jets must be less than 0.85xEb eam . Figure 4a 
shows the distribution of the sum of the energies of hadronic jets. 

After applying the selection, 76 events remain in the data while 73.1 ± 1.4 background events 
are expected. 

For the mode where both W bosons decay hadronically, events satisfying the following 
criteria are selected: 

• The event does not contain any isolated electron or muon and the energy of each non- 
isolated electron or muon must be less than 30 GeV. 

• The number of hadronic jets is at least 4. 

After application of the previous cuts, the dominant remaining background is W-pair pro- 
duction. However, for the signal events, the total visible energy is less than ^/s and the two W 
bosons are not back-to-back due to the energy and momentum carried away by light neutrinos. 
The following additional requirements are then imposed: 

• The visible energy must be less than 0.9x-\/s, as shown in Figure 4b. 

• The acollinearity angle between the two W candidates is greater than 15°, as shown in 
Figure 4c. 

• The total transverse momentum of jets must be greater than 10 GeV, as shown in Figure 
4d. 

For charged heavy lepton masses greater than the W boson mass, the leptons decay to a real 
W. Therefore, for the case where the charged lepton mass is greater than 80 GeV, a kinematic 
fit is applied to improve the determination of jet energies and angles. Both jet-jet invariant 
masses are constrained to the W mass. All possible jet-jet combinations are considered and the 
one which yields the smallest x 2 °f the fit is chosen. 

After applying the selection, 79 events remain in the data while 75.0 ±1.5 background events 
are expected. 



Decay into a stable neutral heavy lepton, L ± — ►L°W ± 

The search for a charged heavy lepton decaying into a stable neutral heavy lepton assumes a 
mass of the associated heavy neutrino L° to be greater than 40 GeV. This assumption is based 
on LEP results at the Z resonance L 3j and implies a large missing energy and large transverse 
momentum imbalance. In the limit of a vanishing mass difference between the charged lepton 
and its associated neutral lepton, Am = m\± — m^o, the signal efficiency is affected by the 
trigger efficiency and the rejection of the two-photon background. Here the search is restricted 
to 5 GeV < Am < 60 GeV. The case of a light neutrino, Am = rriL±, is considered in the 
previous section. The main background is the two-photon process for small mass differences, 
Am < 20 GeVand the qq(7) and W + W~ processes for high mass differences, Am > 20 GeV. 

This signature is very similar to that of a chargino decaying into a stable neutralino and a 
W boson. Therefore, a selection developed for the chargino search [23J is used. It is based on 
the signatures of missing energy, transverse momentum imbalance and isolated leptons. Cuts 
on missing mass and acoplanarity are also used. The selection is optimized for three different 
Am regions: for very low Am region around 5 GeV, the low Am range from 10 to 30 GeV and 
the medium Am range from 40 to 70 GeV. 

After applying the above three selections, 169 events are left in the data while 180.1 ± 6.3 
events are expected from background. 

Search for stable charged heavy leptons 

The analysis follows closely the procedure described in References 0] and E] and extends the 
search to higher masses. The signal events are characterized by two back-to-back high mo- 
mentum tracks in the central tracking detector. Events satisfying the following criteria are 
selected: 

• The event contains two charged tracks, each with momentum greater than 5 GeV and 
polar angle | cos^l < 0.82. 

• The acollinearity angle between the two tracks is less than 15°. 

The momentum and acollinearity angle cuts reduce the background from two-photon pro- 
duced lepton pairs as well as from dilepton annihilation events with a high energy photon in 
the final state. 16,598 events are selected in the data samples with 16,715 events expected from 
Standard Model processes, mainly from e + e~ — > e + e~. 

Stable charged heavy leptons are also expected to be highly ionizing in the tracking chamber. 
The dE/dx measurement is calibrated with Bhabha scattering events, and the mean value of 
the dE/dx distribution is normalized to one. Its resolution is 0.08 units. The ionization energy 
loss has been studied using low energy protons and kaons in our hadronic events up to 10 units 
in dE/dx which corresponds to va\^/E^ earn = 0.995. Figure 5 shows the dE/dx measurements 
of the tracks of data events passing the previous cuts, as well as the simulated signal from pair- 
production of stable charged heavy leptons. We select candidate events for which the ionization 
energy loss for each track is between 1.25 and 8.0 units and the product of the track ionization 
losses is larger than 2.0 units. The upper cut at 8 units rejects very highly ionizing particles 
for which saturation effects may lead to track reconstruction inefficiencies. 

Three candidate events satisfy the selection requirements in the data. The candidate events 
were recorded at beam energies of 97.8 GeV, 102.7 GeVand 103.3 GeVand correspond to pair 



production of stable charged particles of masses 82.9 ±1.8 GeV , 86.5 ±1.8 GeV and 86.2 ± 
2.1 GeV, respectively. However, no candidate events are observed that are consistent with a 
lepton mass larger than 90 GeV. 

The background is estimated using a calibration sample taken on the Z resonance during 
the 1999 and 2000 data taking periods. From the number of two track events from this sample 
passing all stable heavy lepton selection cuts, the background is estimated to be 4.1 ±1.8 events 
for the entire high energy data sample. 

Results 

No signal for new heavy leptons is observed. The number of selected events and expected 
backgrounds are in good agreement. The systematic uncertainty, which is mainly due to the 
uncertainties in the energy calibration, the lepton identification efficiency and purity, and the 
Monte Carlo statistics, is estimated to be 5% relative. To obtain exclusion limits, the selection 
efficiency is reduced by the total systematic uncertainty. These results are combined with our 
earlier L3 data recorded &t y/s = 133 - 189 GeV |U|2l|. 

The selection efficiencies are determined by Monte Carlo simulation. For unstable neutral 
heavy leptons in the mass range from 90 to 100 GeV, the efficiency is 32% to 40% for the 
electron decay mode, and 27% to 35% for the muon decay mode. For the tau decay mode, the 
selection efficiency ranges from 13% for a heavy lepton mass of 80 GeV to 20% for a mass of 95 
GeV. The selection efficiency for unstable charged heavy leptons of mass 100 GeV decaying into 
light neutrinos is 24%. In the case of the charged lepton decaying into its associated neutral 
lepton, the selection efficiency is 5% to 40% with increasing mass difference Am. 

The selection efficiency for stable charged heavy leptons is a function of rn^/E beam and ranges 
from 60% to 70% over the range 0.8 to 0.99, including the trigger efficiency. For m L /E beam > 
0.99, the efficiency drops due to the upper cut on track ionization. 

Taking into account the luminosities, the selection efficiencies and the production cross sec- 
tions we calculate the 95% confidence level, CL, lower limits on the masses of heavy leptons [25J, 
which are shown in Tables 1 and 2. The limits for vector heavy leptons are higher than for 
the sequential and mirror ones due to their production cross section. Figure 6 shows the 95% 
CL exclusion contour in the m L ± — m L o mass plane for decays of the charged lepton into its 
associated neutral lepton. The lower limit on the mass of the charged heavy lepton is shown in 
Table 2 for a mass difference greater than 15 GeV. These limits include the full LEP data set 
and improve upon and supersede our previous results. 
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Table 1: 95% CL lower mass limits in GeV for pair-produced neutral heavy leptons. 



Decay mode 


Model 


Dirac 


Majorana 


L°^eW 


Sequential 
Vector 
Mirror 


101.3 

102.6 
100.8 


89.5 
89.5 


L°^/iW 


Sequential 
Vector 
Mirror 


101.5 
102.7 
101.0 


90.7 
90.7 


L°^rW 


Sequential 
Vector 
Mirror 


90.3 
99.3 
90.3 


80.5 
80.5 



Table 2: 95% CL lower mass limits in GeV for pair-produced charged heavy leptons. 



Decay mode 


Model 




I± ^z/W 


Sequential 
Vector 
Mirror 


100.8 
101.2 
100.5 


I± _► L°W 


Sequential 
Vector 
Mirror 


101.9 
102.1 
101.9 


Stable 


Sequential 
Vector 
Mirror 


102.6 
102.6 
102.6 
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Figure 1: The distribution of the energy in a 30° cone around the second most energetic electron 
candidate. The points are the data, the solid histogram is the background Monte Carlo. The 
shaded histogram represents the simulated signal for e + e~ — > L°L° for a sequential Dirac neutral 
lepton mass of 101 GeV. Both histograms are normalized to the luminosity of the data. The 
arrow indicates the value of the applied cut; all other selection cuts are applied. 
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Figure 2: The distribution of the sum of the muon energies, E^i + E^, normalized to the beam 
enery. The points are the data, the solid histogram is the background Monte Carlo. The shaded 
histogram represents the simulated signal for e + e~ — > L°L° for a sequential Dirac neutral lepton 
mass equal to 101 GeV. Both histograms are normalized to the luminosity of the data. The 
arrow indicates the value of the applied cut; all other selection cuts are applied. 
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Figure 3: The distribution of the cosine of the polar angle of the missing momentum, cos8 m i SS . 
The points are the data, the solid histogram is the background Monte Carlo. The shaded 
histogram represents the simulated signal for e + e~ — > L°L° for a sequential Dirac neutral 
lepton mass equal to 80 GeV. The normalization for the signal Monte Carlo is scaled by a 
factor of 2 for better visibility. The arrows indicate the values of the applied cut. 
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Figure 4: The distribution of a) the sum, Ej etl + Ej et2 , of the energies of the hadronic jets 
normalized to the beam energy, b) the normalized visible energy in the event, c) the acollinearity 
angle of the two W candidates, d) the total transverse momentum, p t . The points are the data, 
the solid histogram is the background Monte Carlo. The shaded histogram represents the 
simulated signal for e + e~ — > L + L _ for a sequential charged lepton mass equal to 95 GeV. The 
normalization for the signal Monte Carlo is scaled by a factor of 2 for better visibility. The 
arrows indicate the values of the applied cuts after all the other selection requirements. 
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Figure 5: The two dimensional distribution of the ionization energy loss of the most energetic 
track, dE/dxi, and the least energetic track, dEjdxi- The solid circles represent the data 
and the open circles represent the simulated signal for e + e~ — » L + L~ for a lepton mass of 
0.985 x Ebeam- The signal normalization is arbitrary. The lines indicate the values of the 
applied cut. 



17 



100- 



<D 80- 



60- 



40 



40 



Sequential and mirror 
Vector 



L3 



L^L°W f 




Excluded 



i i 



i i r 



1 i i 



60 80 

m L ±(GeV) 



100 



Figure 6: The 95% CL limit in the m^ — m^o plane for the case of a charged heavy lepton 
decaying into a stable neutral heavy lepton. The solid and dashed lines represent respectively 
the sequential and mirror, and vector types of heavy lepton. 
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